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Abstract 

A lattice model of critical dense polymers is solved exactly on a cylinder with finite circum- 
ference. The model is the first member CM(1, 2) of the Yang-Baxter integrable series of 
logarithmic minimal models. The cylinder topology allows for non-contractible loops with 
fugacity a that wind around the cylinder or for an arbitrary number £ of defects that propa- 
gate along the full length of the cylinder. Using an enlarged periodic Temperley-Lieb algebra, 
we set up commuting transfer matrices acting on states whose links are considered distinct 
with respect to connectivity around the front or back of the cylinder. These transfer matrices 
satisfy a functional equation in the form of an inversion identity. For even N, this involves a 
non-diagonalizable braid operator J and an involution R = —(J 3 — 12 J)/ 16 = (— 1) F with 
eigenvalues R = (— l)^ 2 . This is reminiscent of supersymmetry with a pair of defects inter- 
preted as a fermion. The number of defects £ thus separates the theory into Ramond (£/2 
even), Neveu-Schwarz (£/2 odd) and Z 4 {£ odd) sectors. For the case of loop fugacity a = 2, 
the inversion identity is solved exactly sector by sector for the eigenvalues in finite geometry. 
The eigenvalues are classified by the physical combinatorics of the patterns of zeros in the 
complex spectral-parameter plane. This yields selection rules for the physically relevant solu- 
tions to the inversion identity. The finite-size corrections are obtained from Euler-Maclaurin 
formulas. In the scaling limit, we obtain the conformal partition functions as sesquilinear 
forms and confirm the central charge c = —2 and conformal weights A, A = A t = (t 2 — l)/8. 
Here t = £/2 and t = 2r — s G N in the £ even sectors with Kac labels r = 1, 2, 3, . . . ; s = 1, 2 
while t G Z — | in the £ odd sectors. Strikingly, the £/2 odd sectors exhibit a W-extended 
symmetry but the £/2 even sectors do not. Moreover, the naive trace summing over all £ 
even sectors does not yield a modular invariant. 
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1 Introduction 



Familiar materials such as plastics, nylon, polyester and plexiglass are made from polymers. 
Polymers [Ij consist of very long chain molecules with a large number of repeating structural 
units called monomers. Polymers exist in low- or high-temperature phases which are char- 
acterised as either dense or dilute. Polymers are dense if they fill a finite (non-zero) fraction 
of the available volume in the thermodynamic limit. 

The modern era of two-dimensional polymer theory began in the late eighties [21 El H] 
when Saleur and Duplantier initiated the study of polymers as a conformal field theory 
(CFT). From the viewpoint of lattice statistical mechanics, polymers are of interest as a 
prototypical example of a system involving (extended) non-local degrees of freedom. It 
might be expected that the non-local nature of these degrees of freedom has a profound 
effect on the associated CFT obtained in the continuum scaling limit. Indeed, the associated 
CFT is in fact logarithmic [5] in the sense that, for certain representations, the Virasoro 
dilatation generator L is non-diagonalizable and exhibits non-trivial 2x2 Jordan blocks. 

In this paper, we consider a Yang-Baxter integrable model of critical dense polymers on 
a cylinder, both on the lattice and in the continuum scaling limit. In fact, this model is the 
first member CM(1, 2) of the infinite series of logarithmic minimal models £Ai(p,p') [6\. An 
alternative approach to logarithmic CFT based on quantum spin chains appears in |7j. The 
fugacity of contractible loops for £A4 (1, 2) is = 0. Previously [8], we considered this model 
on finite-width strips with various boundary conditions. The associated integrals of motion 
and Baxter's Q matrix have been considered in |9J. The solvable critical dense polymer 
model on a cylinder is built from a (locally) planar version [10] of the periodic Temperley- 
Lieb algebra [TTJ. Because of the non-local degrees of freedom, logarithmic theories are 
sensitive to the topology. So changing the topology from a strip to a cylinder has profound 
effects. Specifically, a cylinder topology allows for non-contractible loops with fugacity a that 
wind around the cylinder or for defects that propagate along the full length of the cylinder. 
These can have dramatic effects on the properties of the model. The number of defects £ is 
a quantum number that separates the theory into Ramond (£/2 even), Neveu-Schwarz (£/2 
odd) and Z 4 {£ odd) sectors. Following [U [T2], we use the terminology of supersymmetry 
even though we do not claim any superconformal symmetry in our model. 

Remarkably, as for the square lattice Ising model [13] , the commuting single-row transfer 
matrices of this model satisfy a simple functional equation in the form of an inversion identity. 
This enables us to solve for the exact eigenvalues of the transfer matrices on a finite lattice 
for a = 2. The conformal spectra are readily accessible from finite-size corrections. In 
particular, in the continuum scaling limit, we obtain the partition functions as sesquilinear 
forms and confirm the central charge c = —2 and conformal weights A, A = A t = (t 2 — 1)/8. 
Here t = £/2 and t = 2r — s G N in the £ even sectors with Kac labels r = 1, 2, 3, . . . ; s = 1, 2 
while t e Z-i in the £ odd sectors. On the strip, this model admits a W-extended conformal 
algebra and is identified [H] with symplectic fermions [15] . Strikingly, on the cylinder, we 
find that this extended symmetry only holds in the sum over £/2 odd sectors and not in the 
sum over £/2 even sectors. Moreover, the naive trace summing over all even N sectors does 
not yield the known modular invariant [16] of the c\^ triplet model. 

The layout of this paper is as follows. In Section El we define the solvable critical 
dense polymer lattice model. We also discuss the periodic Temperley-Lieb (TL) algebra, its 
enlargement by adding the shift operators Q and and its relation to the cylinder TL 
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algebra which is a direct generalization of the planar TL algebra of Jones [10] . In Section [31 
we define the single-row transfer matrices directly in the cylinder TL algebra. We also define 
the vector spaces of link states on which these transfer matrices act and relate these to 
the cases of distinct (DC) and identified (IC) connectivities [T7]. The inversion identities 
are derived in Section HJ first in the setting of the cylinder TL algebra and then as matrix 
inversion identities. Certain details are deferred to Appendix |A] The inversion identities 
are solved sector by sector for the transfer matrix eigenvalues in Section [5j while the finite- 
size corrections are extracted in Section [6j In Section [7J the physically relevant solutions 
are obtained empirically and encoded by applying physical combinatorics supplemented by 
selection rules. Finitized conformal partition functions are obtained as sesquilinear forms in 
finitized characters [18J. The conformal partition functions arising in the continuum scaling 
limit are discussed in Section [HJ We conclude with some remarks and directions for future 
research in Section [9j 



2 Lattice Model 

2.1 Critical dense polymers 

To model critical dense polymers on a cylinder with a finite circumference, we consider a 
square lattice on a strip, with N columns and M rows of faces, and identify the left and 
right edges as shown in Figure 1. 




(2.1) 



Figure 1: A typical dense polymer configuration on a 10 x 4 cylindrical lattice. The left 
and right edges are identified to form a cylinder. In this case, the circumference is iV = 10 
and M = 4. No local closed loops are formed on the surface of the cylinder in accord with 
the vanishing loop fugacity = 0. It is, however, possible to allow non-contractible loops 
encircling the cylinder with fugacity a ^ 0. 

An elementary face of the lattice can assume one of two configurations with different 
statistical weights 



or 



(2.2) 



where the arcs represent local segments of polymers. The two possible configurations can be 
combined into a single face operator as 



X(u) 



u 



COSM 



+ sinw 



(2.3) 
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where u is the spectral parameter related to spatial anisotropy. The lower left corner is 
marked to fix the orientation of the square. 

Since the polymer segments pass uniformly through each face, this is a model of dense 
polymers — in the continuum scaling limit, a polymer is space-filling and has fractal dimen- 
sion 2. The non-local degrees of freedom correspond to a number of polymers. It is often 
convenient to think of these degrees of freedom as non-local connectivities. 

Critical dense polymers corresponds to the first member £.M(1,2) of the infinite series 
CAi (p, p') of logarithmic minimal models [6] . Each logarithmic model is characterized by a 
crossing parameter A = n related to the loop fugacity (3 by 

p' — p 

(3 = 2 cos A = 2 cos n, p,p' coprime (2.4) 

pi 

It was argued in [6] that the scaling limits of these integrable lattice models yield logarithmic 
CFTs. In the case of critical dense polymers, 

A = - = (2.5) 

implying that local contractible loops are not allowed. 



2.2 Cylinder Temperley-Lieb algebra 



For the purposes of this paper, the cylinder Temperley-Lieb (TL) algebra is a diagrammatic 
algebra built up from elementary faces. The faces are connected such that the midpoint of 
an outer edge of a face, called a node, can be linked to a node of any other (or even the 
same) face as long as the total set of links make up a non-intersecting web of connections 
on the surface of a cylinder. The cylinder TL algebra is equivalent to the annular algebra of 
Jones [TO] . 

Two basic local properties of the cylinder TL algebra are the inversion relation 




and the Yang-Baxter equation (YBE) [T3 



cos 2 f 



(2.6) 
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u 



u 






{ u 




V 






u — V 



(2.7) 



These are identities for 2- and 3-tangles, respectively, where a fc-tangle is an arrangement 
of faces with 2k free nodes. The identities are established by writing out all the possible 
configurations, while keeping track of the associated weights, and collecting them in classes 
according to connectivities. The left side of (|2.6|) . for example, thus corresponds to a sum of 
four terms of which one vanishes since (3 = 0. The remaining three terms fall into the two 
connectivity classes 





and 



<8s$> 



(2.8) 



The weights accompanying the two equivalent configurations cancel since cosvsin(— v) + 
sin?; cos (— v) = 0, while the last diagram comes with the weight cosf cos(— v) thereby yield- 
ing the identity (12.61) . 

Particular elements of the cylinder TL algebra are the shift or winding operator Q and 
its inverse Q~ l 



n 



-1 



n 



The multiplication implied in 























L 


n- 1 


n = 


i 







(2.9) 
(2.10) 

(2.11) 



is vertical concatenation, and / is the (vertical) identity operator linking every node on the 
upper horizontal edge to the node directly below it on the lower horizontal edge. 



2.3 Periodic Temperley-Lieb algebra and its enlargement 

The periodic (affine) TL algebra [H] of size N is generated by the identity I and the gener- 
ators ej 

TC{N) = (I, ea, e!,..., ex-i) (2.12) 
subject to the periodicity constraints 

ej = e jmodN , j eZ (2.13) 

and the relations 

e) = (3ej, ejej^ej = e jt j = 0, 1, . . . , N - 1 (2.14) 

e^efc = e^ej, j — k ^ 0, ±1 mod iV (2-15) 

where /3 is the fugacity of contractible loops. The TL generators ej are represented diagram- 
matically by monoids (20] . The face operator at position j is given by 

Xj(u) = cos u I + sin u ej (2-16) 

and obtained by rotating the face operator in ( 12. 3ft by 45 degrees in the counter-clockwise 
direction. The periodic TL algebra is infinite dimensional. 
For iV even, one introduces the combinations 

n N rt N 

E = e e 2 e 4 ---e^_ 2 , F = ei e 3 e 5 • • • e N . u E = (3~E, F = (3~F (2.17) 
Letting a denote the fugacity of non- contractible loops, these can be 'removed' in pairs 

EFE = a 2 E, FEF = a 2 F (2.18) 
For iV odd, non-contractible loops cannot appear. 



5 



The elements Q and Q _1 (12.101) of the cylinder TL algebra are not elements of the 
TL algebra. They may be included [19J . though, thereby enlarging the periodic TL 
Translation on the TL generators is then implemented by conjugation 



e 7 _i 



so that the enlarged TL algebra is generated by three independent generators 

STC(N) = (eo^lT 1 ) 

The periodicity constraints ( 12 . 13[) and the relations (I2.14p and (12.151) now read 

Q N e Vl- N = e 
e\ = f3e , eo^eoO^eo = e 
eoQieoQ^i = Q^eoQ^^eo, j = 2, 

For N even, we furthermore have 



iV-2 



e = (e n- 2 fn N = n- N (n 2 e f 

F = fi- 1 (e fi- 2 ) f fi iV+1 = fi- 7V - 1 (fi 2 e ) f fi 
E = ttFtt- 1 , F = tt^Ett 



and 



EVL ±l E = aE, FQ ±X F = aF 
indicating that non-contractible loops can now be removed one by one. 



periodic 
algebra. 



(2.19) 

(2.20) 

(2.21) 
(2.22) 
(2.23) 

(2.24) 

(2.25) 
(2.26) 

(2.27) 



3 Cylinder Transfer Matrix 
3.1 Single-row transfer matrix 

Having introduced the cylinder TL algebra, we now define diagrammatically the single-row 
iV-tangle 



T(u) 



u 














u 
















> 



(3.1) 



consisting of N faces where the dependence on N is suppressed. The left and right edges 
are identified in accord with the periodicity of the cylinder. As discussed below, T(u) has a 
natural matrix representation when acting vertically from below on a given set of periodic 
link states. We thus refer to it as the (single-row) "transfer matrix" , even though it is defined 
as a cylinder iV-tangle without reference to any matrix representation. The shift operator 
Q and its inverse enter naturally as the limits 



limT(w) = tt, lim T(u) = IT 1 

u— s-0 m— >A 



(3.2) 
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Using standard diagrammatic arguments [13J, it follows that T{u) gives rise to a com- 
muting family of transfer matrices where 

[T(u),T(v)} = (3.3) 

As in (12.111) . the implied multiplication in the commutator means vertical concatenation of 
the two iV-tangles in the cylinder TL algebra. It follows, in particular, that [T(u), Q^ 1 ] = 0. 



3.2 Hamiltonian limit 

The Hamiltonian limit of the transfer matrix T(u) is defined in the cylinder TL algebra as 
the iV-tangle appearing as the first sub-leading term in an expansion with respect to u. We 
define H as this iV-tangle up to a factor of Q, that is, 



T(u) = n[l - uH + 0{u A 
where I is the vertical identity diagram 



(3.4) 



(3.5) 



It follows that 
// = 



+ ... + 











\ 






r 



(3.6) 



which in terms of the generators of the 'periodic TL algebra merely corresponds to 

N-l 

H = - ej (3.7) 

j=0 



3.3 Link states 

A matrix representation of T(u) is obtained by acting with T(u) from below on a suitable 
vector space of link states. Suppose there are N nodes arranged periodically around the 
upper horizontal edge of the cylinder. For N even, a link state specifies how these N nodes 
are linked together. Two nodes can be connected by the front of the cylinder or by the 
back. We can consider these two connections as distinct or we can choose to identify the 
two connections and their corresponding link states. In the latter case, we can think of a 
hemi-spherical cap placed on the top of the cylinder so that a connection by the back of 
the cylinder can be continuously deformed to a connection by the front and vice versa. The 
link states are locally planar in the sense that connections are not allowed to cross on the 
extended surface of the cylinder (or capped cylinder). 

The enlarged TL algebra is the appropriate algebra when acting on link states with 
distinct connectivities (DC). In the topology associated with the action on link states with 
identified connectivities (IC), all loops become contractible. This implies that the appropriate 
algebra in this case is the enlarged TL algebra with a — /3. In either case, the enlarged TL 
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algebra is effectively finite when acting on the (DC or IC) link states since the corresponding 
matrix realizations satisfy 

n N = (ir 1 )^ = i (3.8) 

This means that we do not keep track of windings around the cylinder. As matrices, the 
inverse shift operator is given by the Hermitian conjugate of the operator itself 

Q- 1 = Q j (3.9) 

For N = 4, there are six DC link states 

f^r^\ , \ r\ r , ^\ fr , r\ r\ , *\ fr\ , C (3.10) 

and two IC link states 

(^\ , rs (3.11) 

In general, the dimension of the vector space of link states for DC and IC is given by the 
central binomial coefficients and Catalan numbers 



dim(V^) = m, dim(y^) = — I j, n = -, N even (3.12) 

A node that is not linked to another node gives rise to a defect which may be viewed as a 
link to the point (above) at infinity. For A^ odd, there is at least one defect. In the presence 
of defects, there is only one way (either by the front or by the back) to connect two nodes, 
so there is no distinction between identified and distinct connectivities. The dimension of 
the space of link states with precisely I defects is 

dim(V^ ) ) = UjJS ' £ = N mod 2 ( 3 - 13 ) 

It is the set of DC link states which corresponds to this for £ = 0: = V^ c . For A^ = 3 
and i — 1, there are 3 link states 

I rs , rs I , -\ I r (3.14) 

We note that (13. 8p and (13.91) remain valid when acting on link states with defects. Also, 
defects can be annihilated in pairs, but not created, by the action of the cylinder TL algebra. 
If we allow for an arbitrary number of defects of given parity, the number of link states is 

=*"-'. E(£l)^--i(f) (3.15) 
3.4 Augmented link states 

When the direction of transfer is fixed, the cylinder TL algebra reduces to the enlarged 
(periodic) TL algebra. The transfer A^-tangle then acts naturally on link states with Af 
nodes. For the purpose of computer calculations, however, the transfer matrix is conveniently 
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written in terms of the enlarged TL algebra £TC(N + 2) acting on a suitable vector space 
of link states of size N + 2. Explicitly, 



N-2 



T(u) 



e N 



I I ^(u)n 



(3.16) 



as shown in Figure [2], where the face operators are given by f)2.16p . The link states are all 
augmented by a spectator half-arc joining the nodes in positions iV and N + 1. For N = 4, 




N N+l 



Figure 2: Diagrammatic representation of the single-row transfer matrix built from the 
enlarged perodic TL algebra acting on iV + 2 strings. 



for example, the six augmented DC link states are given by 



(3.17) 



where the nodes are labelled from j = to j = N + 1. For N = 3, and an arbitrary odd 
number of defects, we have four augmented link states 



1^^ , , I (r\ , 111^ (3.18) 

where the nodes are labelled, as before, from j = to j = N + 1. 

To further facilitate the implementation of the analysis on a computer, we note that the 
link states can be described as sets of pairs of connecting arcs and defects. In this language, 
the set of link states f)3.14p reads 

{{1},{2,3}},{{1,2},{3}},{{2},{3,1}} (3.19) 



9 



4 Inversion Identities 



Remarkably, the transfer matrix (13. ip satisfies an inversion identity in the cylinder TL alge- 
bra. This identity is thus independent of the choice of vector space of link states eventually 
acted on to form a matrix representation. First, we describe the identity in the general 
cylinder setting (and prove it in Appendix lA.lj) . We then characterize it when acting on 
the various link states introduced above. The inversion identity is unique to critical dense 
polymers among the logarithmic minimal models [;6j. The analogous inversion identity for 
critical dense polymers on the strip is discussed in [8]. Although the inversion identities on 
the strip and the cylinder have common features, their solutions and general properties are 
very different. 



4.1 Cylinder inversion identity 

In preparation for the inversion identity, we introduce the two 3-tangles 

— -f- and X- 



(4.1) 



as they play important roles as building blocks in the following. The Af-tangle «7, in particu- 
lar, is defined as the sum of the 2 N possible horizontal combinations of A" of these 3-tangles, 
where the left and right edges are identified to respect the cylinder topology. Due to the 
minus sign in (14. ip . exactly half of the terms in J appear with a minus sign. For small N, 
we thus have 



|JV=1 



\N=2 



\N=3 



2 + £ 



-ft 2 + ft" 



-ft 2 



ft- 



(4.2) 
(4.3) 



m 

m-m-m 



ft 



-2 



(4.4) 



Although the A^-tangle J = J (a) depends on the fugacity a of non-contractible loops, it is 
independent of the spectral parameter u. 



Inversion Identity The A '-tangle T(u) defined in ( 13. lj) satisfies 



7T , 



T(u)T(u+-) 



(cos 2N u + (-l) N sm 2N u)l+{cosusmu) N J (4.5) 

The proof of this Inversion Identity is provided in Appendix IA.1I 

It follows from the inversion identity (14. 5 p and the commutativity property (13. 3 p that J 
is a symmetry of the model in the sense that it commutes with the transfer matrix, 

[J,T(u)] = [J,^ 1 ] = (4.6) 

As discussed in Section 14. 3[ J is related to the so-called braid transfer matrices. 
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4.2 Matrix realization of J 



Obtaining a matrix realization of J is greatly simplified by the Drop-Down Lemma below. 
To state it, we introduce the arc-part of a link state as the part remaining when ignoring 
all defects. For given number of nodes N, a (DC or IC) link state is thus characterized 
completely by its arc-part and its number of defects. The arc-part of any of the three link 
states in ( 13. 14j) consists of a single half-arc whose position depends on the original link state. 
We also say that the arc-part of a link state is contained in the arc-part of another link state 
if the bigger arc-part can be constructed from the smaller one by addition of half-arcs. Two 
identical arc-parts are said to be contained in each other. 

Drop-Down Lemma The action of J on a given input link state results in link states 
whose arc-part contains the arc-part of the input link state. 

The proof of this Drop-Down Lemma is provided in Appendix IA.2I 

So far, we have not specified the class of link states which J is acting on. The following 
Sector Lemma concerns the matrix realization of J in a given sector. 

Sector Lemma In a given sector defined by a specified number of defects £, the matrix 
realization of J is diagonal and given by 



The proof of this Sector Lemma is provided in Appendix IA.31 

It is also of interest to examine the action of J on the set of link states with an arbitrary 
number of defects. After completion of the drop-down process, the remaining part of the 
link state consists of defects only and the situation is equivalent to a scenario with system 
size Ni — £ where £ is the number of defects of the original input link state. That is, all the 
essential data is encoded in the Drop-Down Lemma and the action of J on link states with 
defects only. It is noted that this is true for all sectors or combinations thereof, in particular 
for the union of all sectors of the parity of N. 

As discussed in Appendix \AA\ the matrix realization of J acting on the set of link states 
with an arbitrary odd number of defects is the zero matrix. For N even, on the other hand, 
the matrix realization of J acting on the set of link states with an arbitrary even number 
of defects is not diagonal, not even diagonalizable. With respect to the number of defects, 
the matrix is upper block triangular. The blocks on the diagonal are the same as the ones 
obtained by the sector-by-sector analysis above, while the entries outside these blocks give 
rise to a non-trivial Jordan decomposition. In particular, for DC link states with a = 2 
and N = 8, 10, 12, 14, we observe that Jordan blocks of rank 2 appear, but not of higher 
rank, while there are no rank-2 blocks for the smallest system sizes N = 2,4,6. Assuming 
that no Jordan blocks of rank 3 or higher occur for iV > 16, we conjecture that the minimal 
polynomial identity satisfied by J, valid for all even N, is 




N, £ even 



(4.7) 



(j 2 -4/) 2 = 0, a = 2 



(4.8) 




(4.9) 
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The eigenvalues are R — |J = ±1. We will comment on this involution in Section [91 and 
refer to Appendix IA.4I for additional details on the matrix realization of J when acting on 
DC link states with an arbitrary number of defects. 



4.3 Braid limits 

Let us define the braid limits by 
X(u) 



u± 



lim 



e T 4 



I 7TI 



B 



± 



lim 



T(u] 



u->±ioo sin Ar (w + ~ 



V\ (4-10) 



By taking the braid limit of ( 14. 5p . the matrix J is seen to be simply related to the braid 
transfer matrices B ± 

(B±) 2 = 2I+(±i) N J (4.11) 
Using that J = for A" odd, we find 



±\2 



21 +(-!)- J, N even 
2/, JV odd 



-l)t((B ± ) 2 -21), N even 

odd 



(4.12) 



Assuming the conjectured minimal polynomial identity ( 14. 8 j) for J implies that B = B 1 * 1 
satisfies 

B 4 (B 2 -4/) 2 = 0, AT even; B 2 = 2/, N odd (4.13) 

The eigenvalues of the braid matrices are thus of the form 



S7T 

B = 2 cos— , s = 0,l,2,3,4 



where 



B = 0, ±2, N even; 



B = ±V2, A^odd 



(4.14) 
(4.15) 



4.4 Matrix inversion identities 

Once a matrix representation of T(u) has been fixed, the inversion identity ( 14. 5ft translates 
into a matrix inversion identity. Here, we consider the link states counted in ( 13. 13f) . including 
the DC link states for £ = but not the IC link states. The latter are discussed in Section 1431 
It follows from the Inversion Identity (I4.5p and the Sector Lemma (14. 7p that the matrix 
inversion identity for a given sector reads 

T{u)T{u+^) = (cos 2Af M + sin 2iV M + (-l)^(2 + (a 2 -4)^o)(cosMsin M f)/ (4.16) 

for A^ even, while for A^ odd, it reads 

(4.17) 



T(u)T(u + -) = (cos 2N u-sm 2N u)l 



In order to solve the inversion identity ( 14.16)) for the associated eigenvalues (N even), we 
choose to focus on particular values for a in the following. For a 2 = 4, we thus have the 
matrix inversion identities 



T(u)T(u + ^) = (cos N u+(-l) 



2 sin u) 1 



(4.18) 
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while, for a = and i = 0, we have 

T(«)T(w + |) = (cos^u- (-l^sin^Vfl (4.19) 

We find that 

a = 2 (4.20) 
is the most natural value for the fugacity of non-contractible loops. 

4.5 Identified connectivities 

The distinction between DC and IC link states is only meaningful for link states without 
defects. In this case, the number of IC link states is smaller than the number of DC link 
states, cf. (I3.12p . implying that, for given N, the corresponding matrix realization of T(u) 
or J is of lower dimension in the IC case than in the DC case. As already mentioned, the 
appropriate algebra in the IC case is the enlarged TL algebra with a = j3 = 0. It follows 
that the matrix inversion identity takes the same form as (I4.19p . namely 

T(u)T(u + -) = (cos N u-{-l)%sin N ufl (4.21) 

but applies to matrices of smaller dimension than the ones appearing in (I4.19p . 

5 Solution on a Finite Cylinder 

In this section, we solve the inversion identities for the transfer matrix eigenvalues on finite- 
size cylinders sector by sector for a = 2. We also discuss the relation, through finite-size 
corrections, to the conformal partition functions. The methods build on the previous works 
[H], [21] and [8]. 

The key idea is that the eigenvalues T{u) of the transfer matrices in a given sector are 
determined, up to an overall constant p, by the positions Uj of their zeros in the analyticity 
strip -7r/4 < Reu < 37r/4 

N 

T{u) = p Y\ sin(w — uj) (5.1) 
i=i 

These eigenvalues are Laurent polynomials in z = e m . Given a (right) eigenvector indepen- 
dent of u, this follows since each entry of the transfer matrix is of this form. This argument, 
which applies to the action of the transfer matrix on a particular eigenvector, holds even 
if the commuting transfer matrices are not diagonalizable. Typically, in the I even sectors 
on the cylinder, we find that the transfer matrices are not diagonalizable as we find Jordan 
blocks of rank 2. We have not found Jordan blocks of any higher rank. Moreover, we find, 
in a given I sector, that all of the eigenvectors (excluding generalized eigenvectors) are inde- 
pendent of u. Degenerate eigenvalues are exactly degenerate as Laurent polynomials so such 
eigenvalues occur with degeneracy 2. These observations enable us to numerically obtain all 
of the eigenvalues as Laurent polynomials in z = e m . 

We emphasize that this situation on the cylinder is in contrast to critical dense polymers 
on the strip in sectors of the extended Kac table for which the transfer matrices are found 
empirically [S] to be simultaneously diagonalizable. It seems that it is not possible to avoid 
reducible yet indecomposable representations of rank 2 in the £ even sectors on the cylinder. 
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5.1 Finite-size corrections 



The partition function of critical dense polymers on a periodic lattice of N columns and M 
rows is defined by 



TyT(u 



,A7 



£r„ 

n>0 



U 



n>0 



-M£ n (u) 



(5.2) 



Here the sum is over all eigenvalues of T(u), including possible multiplicities, and S n (u) with 
n = 0, 1,2, . . . is the energy associated to the eigenvalue T n {u). The maximal eigenvalue 
T Q (u) is labelled by n = 0. The maximal eigenvalue in the sector with I defects is denoted by 
Tq^{u). Conformal invariance of the model in the continuum scaling limit dictates 
that the leading finite-size corrections for large N are of the form 



So 



N fbulk 

2m 



ire 



sin^ 



N 
2tt 



6N 
[(A + k)e~^ 



(A + k)e 



id] 



(5.3) 



— [(A + A + k + k) sin^ + i(A - A + k - k) cos^] 



Here fbulk is the bulk free energy per face [8] 



fi 



bulk 



ilog2 



■k/2 



log 



sint 



+ sin 



2u)dt 



(5.4) 



and $ = 2u is the anisotropy angle. The conformal spectrum is determined by the central 
charge c = —2, the conformal weights A, A and the excitations or descendants labelled by 
the non-negative integers k, k. The conformal weights are given by 



t 2 -l 



t E \7L 



A = A t 

where t can be integer or half-integer 

A r>s = A 2r _ s e {-± 0,|,1, f ,...}, reN,s = l,2; iV even 



(5.5) 



A 



A t e{- 



8' 8 

45 77 
32 ' 32 ' 32 ' 32 ' 32 ' 32 



3 5 2 i r , 77 :17 1 



t e Z- f; 



iV odd 



(5.6) 



Here r, s are the Kac labels [5] and t = £/2 where £ is the number of defects. The Kac table 
is shown in Figure |3j 

In the scaling limit, the conformal partition functions are sesquilinear forms in characters 



Z(q) = ^A/^aXa^Xa® 



(5.7) 



A, A 



where the characters are of the form 



A+k 



k=0 



ch rjS (g) 
ch t (q) = 



-c/24+A r , s Q _ 



q' 



nr =1 (i-^ 

-c/24+At 



nr=i(i-^ 
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(5. 
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Figure 3: Kac table of the I even sectors of critical dense polymers. The relevant rows, 
s = 1 and s = 2, label the Neveu-Schwarz and Ramond sectors, respectively. The (r, s) 
representations indicated with a red quadrant are irreducible representations with characters 

XrM = Ch r,s((l)- 



and d&(k) are the degeneracies at level k. The modular nome is 

M 



q = exp(27nr) 



N 
M 



exp[2(7T — 



-Se 



-2iu 



q = exp(— 2-nir), t = exp[— i(n — ■&)] = —5e 2m 



\Q\ 



qq 



exp(— 47r5sin 2u) 



(5.9) 
(5.10) 
(5.11) 



where 5 = M/N is the aspect ratio and Imr > in the physical strip < Rew < |. 
The characters ch r ^ s (q) with s = 1,2 in the N even sectors are the characters of irreducible 
Kac representations [6 J whereas ch t (q) in the iV odd sectors are the characters of generic 
(irreducible) Virasoro modules. 



5.2 Z 4 sectors (N odd, £ odd) 

In the Z4 sectors, the sector-by-sector inversion identity for the eigenvalues is 



7T , 



'/'(// )'/'(// + -) = cos 2N u - sin 2N u 



Factorizing the right side gives 



2N ■ 2N 

cos u — sin u 



-2Niu 



N 



2 2N- 



AN 



(5.12) 



(5.13) 
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Sharing out the zeros to solve the functional equation, gives 

-Niu N 



T(u) 



i) N ' 2 t 



2^-1/2 



n 



e- + ^tan (2j - 1)7r 



AN 



where e 2 



e? = 1. The ordinates of the locations of zeros are 



(5.14) 



Vj 



, iuyu.il — 



— h log tan 



|)tT 



A typical pattern of zeros is shown in Figure HJ 
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(5.15) 



Figure 4: A typical pattern of zeros in the complex w-plane for the Z 4 sectors (N odd, 
i odd). Here, N = 11 and £ = 7. The ordinates of the locations of the zeros itj are 
y.,- = — | log tan ^Tjy ) j = 1,2, iV. At each position j, there is either a 1-string with 
Reitj = 7r/4 or a 2-string with Keuj = — 7r/4, 37r/4. 

We see that these solutions (eigenvalues) satisfy the crossing symmetry 

m-u) = t(u) 



(5.16) 



The choice ej = —1 for a particular j corresponds to an elementary excitation. In principle, 
up to the overall choice of sign e, there are 2^ possible eigenvalues allowing for all excitations. 
However, only (n-A of these solutions actually occur as eigenvalues and these are determined 

V 2 ' 

by selection rules as explained in Section [0 For £ = 1, the largest eigenvalue T i(m) occurs 
for Ej = 1 for all j = 1,2, . . . , iV , that is, there are 2-strings at each position j and no 
1-strings. The patterns of zeros of T(u) are conveniently encoded by introducing pairs of 
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single-column diagrams as shown in Figure [5j The right column corresponds to the 1-strings 
in the lower-half w-plane (associated with q). The left column corresponds to the 1-strings 
in the upper half-plane (associated with q), including the real axis, but rotated through 180 
degrees. Positions occupied by a 1-string are indicated by a solid circle and unoccupied 
positions are indicated by an open circle. For N odd and £ odd, the patterns of zeros of 
T ^(u) are shown in Figure El 
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Figure 5: The patterns of zeros of T(u) in the Z 4 sectors are encoded by introducing pairs 
of single-column diagrams. The right column corresponds to the 1-strings in the lower-half 
u-plane (associated with q). The left column corresponds to the 1-strings in the upper half- 
plane (associated with q), including the real axis, but rotated through 180 degrees. Positions 
occupied by a 1-string are indicated by a solid circle and unoccupied positions are indicated 
by an open circle. The 1-string energies are given by Ej = l(j — |). Here, N — 11, a — 2, 
a = -2 and 1 = 2(a + a) + 1 = 1. 



(iV+l)/2 


o 




(N-l)/2 


o 


o 




o 


o 


3=3 


o 


o 


3 = 2 


o 


o 


3 = 1 


o 






1 = 1 £ = 3 £ = 5 

a = a = a = a = —1 a = a = 1 o 



o 
o 

11 



-2 a = a = 2 a = a 



Figure 6: Groundstate configurations in the Z 4 sectors with £ defects (N, £ odd). The 
quantum numbers a, a are as shown and the conformal weights are (A, A) = (A^/ 2) A«/ 2 ). 

Next, we separate the contribution from zeros in the upper and lower half-planes by 
keeping ej for j = 1, 2, . . . , (N + l)/2 and setting 

= CAf+i-i = ±1, J = l,2,...,^i (5.17) 
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By convention, we treat the zeros on the real axis labelled by j = (N + l)/2 as if they are 
in the upper half-plane. Setting e = ji nj=i gi ves 

N+l JV-1 

T(u) - (e^tan^) ft (*,«- + <«* (5,8) 

3=1 3=1 

We fix n — +1 to ensure that T 0) i(0) = 1 consistent with ( 13. 2p . Taking the ratio of (I5.18P 
with precisely one €j = —1 or e 3 ■, — — 1 to (I5.18P with all €j — e,- — +1, and then taking the 
limit with a fixed aspect ratio 5 = M/N, gives 

/2m _ - f (gf-pTr x M 

hm (- = exp[-(j-i)7rz5e- 2i "] = g ^ (5.19) 

J™..( 1. ) = exp[ (j -I),ae-| = (5.20) 

where 



M,AT->oo \ e 2iu + ^ CQt 



^ = |(j - |), iV, £odd (5.21) 
It follows that the conformal partition function in the Z4 sectors with I defects is 

Ze(q) = {qq)-^ + ^Y,^ J1)l25{t3 ^ ( 5 - 22 ) 

where the sum and the allowed values of ej and ej are determined by selection rules as 
explained in Section [71 Here S(j,k) = Sj^ is the Kronecker delta. The prefactor involving 
the central charge and conformal weights comes from the largest eigenvalue T g(u) with 
1=1 and is obtained by applying Euler-Maclaurin as explained in Section [6j 

5.3 Ramond and Neveu-Schwarz sectors (N even, £ even) 

In the I even sectors, the sector-by-sector inversion identity for the eigenvalues is 

T(u)T(u+^) = (cos N u + (-l^ N -^ 2 sm N u) 2 (5.23) 



Using the identities 

cos N u + (-l) N / 2 sm N u = ^ (e 4 *" + tan 2 ^w* )' i/2 ^ {5M) 

3=1 

^(2-^m^ ^ 



COS^ 



M -(-l)^sin^ = -^zr- II ( e4m + tan2 ^)' £ / 2odd ( 5 - 25 ) 

3=1 
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we see that for the two parities of £/2 



N/2 



(cos^ + (-l)^sin^) 2 = ^n(e 2 -+z6,t a n^^)(e 2 --z ej tan (2j ^ 



3=1 



2N 



N/2 

X 

3=1 



n («~ + i W tan ^^)(e 2 '» - ift tan <B^) (5 .26) 

3=1 

(cos^-(-l)^sin^) 2 = A 2 6 2jV _ 2 J] (e- + ^tan^)(e--^tan^) 

3=1 

V/2-1 

J] (e 2 - + ^/i J tan^)(e 2 --z / i,tan^) (5.27) 



iV/2-1 

X 



where ej = fij = 1. Sharing out the zeros to solve the functional equation, gives 



T(u) 



e(—i)~e~ u -tt / 2iu ■ (2j-l)7r\/ 2m . (2j-l)n 
i=i 

( } 2N1 II + * tan f ) + ^ tan , f odd 



even 
(5.28) 



We see that these solutions (eigenvalues) satisfy the crossing symmetry 




T(§ -u) = T(u) (5.29) 
The ordinates of the locations of zeros are 

^r^, ^/2even, j = 1, 2, . . . , N/2 

yj = { ^ N ( 5 - 3 °) 

e/2 odd, j = l,2,...,iV/2-l 

A typical pattern of zeros is shown in Figure [71 

The choice 6j = —1 or fij = —1 for a particular j corresponds to an elementary excita- 
tion. In principle, up to the overall choice of sign e, there are either 2 N or 2 N ~ 2 possible 
eigenvalues allowing for all excitations. However, only ( n-A of these solutions actually occur 

2 

as eigenvalues and these are determined by selection rules as explained in Section [71 For 
£ = 0, the largest eigenvalue T 0j o(u) occurs for €j = fij = 1 for all j = 1, 2, . . . , N, that is, 
there are 2-strings at each position j and no 1-strings. The patterns of zeros of T(u) are 
conveniently encoded by introducing pairs of double-column diagrams as shown in Figures [8] 
and [TD1 The right column corresponds to the 1-strings in the lower-half w-plane (associated 
with q). The left column corresponds to the 1-strings in the upper half-plane (associated 
with q), including the real axis, but rotated through 180 degrees. Positions occupied by a 
1-string are indicated by a solid circle and unoccupied positions are indicated by an open 
circle. For N even and i/2 even or odd, the patterns of zeros of T ^(u) are shown in Figures [HI 
andHU 
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Figure 7: A typical pattern of zeros in the complex w-plane for the £ even sectors. Here, 
N = 12 and £ = 2. The ordinates of the locations of the zeros Uj are y~j = — \ log tan = 
1,2, ... , N/2 — 1. At each position j, there is either two 1-strings with Reuj = 7r/4, two 
2-strings with real parts KeUj = — vr/4, 37r/4 or one 1-string and one 2-string. A double zero 
is indicated by a black circle, a single zero by a grey circle and an unoccupied position by 
an open circle. 



Next, we separate the contribution from zeros in the upper and lower half-planes. For 
£/2 even, we keep €j and fij for j — 1,2, ... , [(N + 2)/4j and set 

6j = = ±1, j = 1, 2, . . . , LiV/4j (5.31) 

fij = VN/2+i-j = ±1, 3 = 1,2,..., [N/A\ (5.32) 

For £/2 odd, we keep ej and fj,j for j = 1, 2, . . . , [N/4\ and set 

tj = e N/2 . 3 = ±1, j = l,2,..., [(N - 2)/4j (5.33) 
p.j = Hn/2-3 = ±1, i = 1,2,..., L(iV - 2)/4j (5.34) 

By convention, we treat the zeros on the real axis (labelled by j = (N + 2)/4 in the case £/2 
even and N = 2 mod 4 and by j = N/4 in the case £/2 odd and N = mod 4) as if they 
are in the upper half-plane. For £/2 even, setting e = // n}={ 4 ' ^ A? gives 

' -7V» L(iV+2)/4j 
,' 1 

L7V/4J 



v i 2 N ~ 1 V 3 2N J\ J 2N , 

IJV/4J 

(e 2 - + ^ tan (V» + z cot (5.35) 
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Figure 8: The patterns of zeros of T{u) in the Ramond sectors (N even, £/2 even) are encoded 
by introducing pairs of double-column diagrams. The right double-column corresponds to the 
1-strings in the lower-half w-plane (associated with q). The left double-column corresponds 
to the 1-strings in the upper half-plane (associated with q), including the real axis, but 
rotated through 180 degrees. At each position in a double-column diagram, there are 0, 1 or 
2 1-strings. Positions occupied by a 1-string are indicated by a solid circle and unoccupied 
positions are indicated by an open circle. The 1-string energies are given by Ej = (j — |). 
There are L^J^J positions on the left and [^J on the right. Here, N = 22, a = 2, a = 
and^ = 2(a + a) =4. 
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Figure 9: Groundstate configurations in the Ramond sectors with I defects (N, £/2 even). 
Here N = 14, the quantum numbers a, a are as shown and the conformal weights are 
(A, A) = (A i/2 ,Ae /2 ). 
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Figure 10: The patterns of zeros of T(u) in the Neveu-Schwarz sectors (N even, £/2 odd) 
are encoded by introducing pairs of double-column diagrams. The right double-column 
corresponds to the 1-strings in the lower-half w-plane (associated with q). The left double- 
column corresponds to the 1-strings in the upper half-plane (associated with q), including 
the real axis, but rotated through 180 degrees. At each position in a double-column diagram, 
there are 0, 1 or 2 1-strings. Positions occupied by a 1-string are indicated by a solid circle 
and unoccupied positions are indicated by an open circle. The 1-string energies are given by 
Ej = j. There are [^J positions on the left and L^r^J on * ne right- Here, N = 24, a = 2, 
a = and t = 2(<r + a + 1) = 6. 
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Figure 11: Groundstate configurations in the Neveu-Schwarz sectors with £ defects (N even, 
£/2 odd). Here iV = 16, the quantum numbers a, a are as shown and the conformal weights 
are (A, A) = (A e/2 ,A e/2 ). 
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For i/2 odd, setting e = /i Yij=i e i fa gi yes 

T(«) = M Z) II ( e 2iU + -itana(e ie - + zcot^ 

i=i 



L(iV-2)/4j 

X 

3=1 



n { e2iu + ^ tan ^) (^ e2m + * cot ^) ( 5 - 36 ) 



Fixing n = +1 ensures that T (0) = T 0j2 (0) = 1 consistent with (13. 2p . For £/2 even, taking 
the ratio of (I5.35P with precisely one ej,ej or fij, fij = — 1 to (I5.35P with all €j = fij = e,- = 
ftj = +1 and taking the limit M, N — > oo with a fixed aspect ratio 5 = M/N gives 

/ p 2iu _ - t (2j-l)7r x M 

hm (- 2 ^ = exp[-(2j-l)vr^e- 24 "] = ^ (5.37) 

_ e 2m + icot (2tl)ZL\^ 

2^ 

2./V 



lim (— ' 2 ff = exp[(2j-lWe 2m ] = f (5.38) 



where Ej = (j — |). Similarly, for £/2 odd, taking the ratio of (I5.36P with precisely one 6j 
or fij = — 1 to (15.361) with all 6j = fij = Ej = fij = +1 and taking the limit M, N — > oo with 
a fixed aspect ratio 5 = M/N gives 

l im * tan,JV ) = exp[-2j7ri5e- 2iu ] = (5.39) 



-e + z cot 



lim ^ = exp[2j7rz(5e 2w l = f (5.40) 

M,iv^oo V e 2i « + % cot f 1 1 J 



where E'j = j. 

We conclude that, in the £ even sectors, 



£. = )3~h N even, £/2 even 

S j, N even, ^/2 odd 1 ' j 



and the conformal partition functions are 



f i N + 2 1 i N i 



, g even 



(w) _c/24+Al Xl^ E ' i=i ^ =1 l^i - | odd 



(5-42) 

where the sums and the allowed values of determined by selection rules as 

explained in Section [7J The prefactor involving the central charge and conformal weights 
comes from the largest eigenvalue T ^(u) with i = 0, 2 and is obtained by applying Euler- 
Maclaurin as explained in Section O 
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6 Finite-Size Corrections from Euler-Maclaurin 

In this section, we use the Euler-Maclaurin formula to obtain the finite-size corrections for 
the following three groundstates 

U ~ i ) 7r 

Tofl(u): Ramond (N even, £ = 0) yj = — | log tan — - — double 

Tq^{u): Neveu-Schwarz (N even, £ = 2) yj = —| log tan— double (6.1) 

r ,i(u): Z 4 (N odd, I = 1) y d = -\ log tan 2U ^ ! single 

In each of these cases, there are no 1-strings, only single or double 2-strings at the positions 
given by yj. The maximum eigenvalues in these sectors take the following real forms involving 
only the geometric factor sin 2u 

N N 



T.,.M= ^n(^" + "- 1 ^) 2 = in(-4riE + -2 1I ) (6.2) 

3=1 Z 3=1 bU1 N 

N-2 JV-2 
ATj^-+l p (2-N)iU 2 • 2 yy a -1 

*W«> = II («~ + !tol 4) = £ II (-W + <") 

i=i j'=i # 

• 3iV + 2 —]\[iu N (n • -|\ 1 I 

3=1 3=1 2AT 

It is confirmed that these eigenvalues satisfy T ^(0) = 1 for £ = 0, 1, 2. 

The logarithms of the eigenvalues Tq^{u) with £ = 0,1,2 involve sums of terms which 
are singular at both endpoints in the limit iV — > oo. To remedy this, we introduce [21] the 
function 

r 1 i 1 

F(t) = log t(n-t)( + sin 2m) = log* + log(vr - t) + log ( + sin 2m) (6.5) 

L v sint J sint 

with 

lim F(t) = lim F(t) = logvr, - lim F'(t) = lim F'(t) = --sin2u (6.6) 

For simplicity, we suppress the u dependence. The endpoint or midpoint Euler-Maclaurin 
formula [21] can now be applied to approximate the sum by an integral 

M 1 r b i h 

5>(a + *Zi) = r / F(t)rft+-[F(6) + F(a)] + -[F'(6)-F'(a)] + 0(/ i 2 ) (6.7) 

M-l 6 

g F(a + (k + \)h) = - ^ F(t)dt - -[F'(b) - F'(a)] + 0(h 2 ) (6.8) 

where b = a + Mh and h is small. Due to (I6.5p . we also need the asymptotic expansion of 
the logarithm of the gamma function 

logrfe) = (y-|)logy-y + |log(27r) + ^- + 0(r 2 ), y^oo (6.9) 
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First, let us consider the Neveu-Schwarz sector with N even and £ — 2. Setting tj 
2j7r/N, we have 



JV N 
2 2 



-logT 0)2 (u) = f \og2- log N - F(t j) + 2 J2 l °Stj (6-10) 

3=1 3=1 

The sum over F can be approximated using the endpoint Euler-Maclaurin formula, with 
a = 0, b = 7r, h = 2ir/N, and the sum over logarithms by the asymptotics of the gamma 
function. This yields finite-size corrections of the form (15. 3 p with c = —2 and A = A = 
Ax = 0. 

Next, consider the Ramond sector with N = mod 4 and I = 0. Setting tj = (2j — 1)tt / N , 
we have 

V_1 _V_ 

2 2 

- logT ,oH = (f - 1) log2 - F (^) + 2 Y. l °Z t i ( 6 ' n ) 

The sum over F can be approximated using the midpoint Euler-Maclaurin formula, with 
a = 0, b = 7i, h = 2tc/N, and the sum over logarithms by the asymptotics of the gamma 
function. This yields finite-size corrections of the form (15. 3 p with c = —2 and A = A = 
A = -1/8. 

Lastly, consider the Z 4 sector with N odd and I — 1. Setting tj = (2j — 1)tt/2N, we have 



Ni N 
2 2 



-logT » = ^\og2-lJ2 F ^ + J2 l0 ^ ( 6 - 12 ) 

fe=0 j=l 

The sum over F can be approximated using the midpoint Euler-Maclaurin formula, with 
a = 0, b = tt, h = tt/N, and the sum over logarithms by the asymptotics of the gamma 
function. This yields finite-size corrections of the form (15. 3 p with c = —2 and A = A = 
Ai/a = -3/32. " 

In summary, the conformal finite-size predictions are of the expected form (15. 3 p with 

c=-2, A = A = A //a = ~ ~ 0, £ = 0,1,2 (6.13) 

The Euler-Maclaurin analysis can be extended to the excitations in these sectors as in [8] 
but we do not do this here. Instead, we use physical combinatorics to directly obtain the 
energy levels of the excited states. 

7 Physical Combinatorics and Selection Rules 

7.1 Z 4 sectors (JV odd, I odd) 

The building blocks of the spectra in the Z4 sectors consist of the g-binomials 



n 




n 


m 




_[n/2\ - a_ 



i^+i) J2 qZjmsE^ a =\n/2\-m (7.1) 



*? cr-singlc 
columns 
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The sum is over all single- column diagrams as in Figure [TBI with a fixed 



with Ej = ! ( j — 

a. Here a is a quantum number given by the number of 1-strings at even positions j minus 
the number of 1-strings at odd positions j 



a = m evcn - m odd 



L«/2j L("+1)/2J 

^ m 2 k - ^ m2fc -! 

k=l k=l 



The number of 1-strings plus 2-strings at any given position is exactly one 

1, 3 = 1,2, 



Tin 



n 



(7.2) 



(7.3) 



The single-columns with quantum number a are generated combinatorially by starting 
with the minimum energy configuration of 1-strings for given a as shown in Figure [T21 
Empirically determined selection rules dictate that in a sector with £ defects the quantum 
numbers of the groundstate satisfy 



a 



a 



{£ - l)/4, £ = 1 mod 4 
-(£+l)/4, £ = 3 mod 4 



The energy of these groundstates is E(a) + E(a) = jq(£ 2 



14(7+ II 



1). 



1,3,5,7,... (7.4) 




a 







Figure 12: Z 4 sectors (N, £ odd): Minimal configurations of single-columns with energy 
E(a) = \o(o + |). The quantum number cr = \n/2\ — m is given by the excess of blue (even 
j) over red (odd j) 1-strings. At each empty position j, there is a 2-string. This analyticity 
strip is in the upper-half complex w-plane rotated by 180 degrees so that position j = 1 
(furthest from the real axis) is at the bottom. 

Excitations, incrementing the energy by one unit, are generated either by inserting a 
pair of 1-strings at positions j = 1 and j = 2 or incrementing the position j of a 1-string 
by 2 units. The g-binomials are illustrated in Figure [T3J Empirically, we find that all the 
excitations satisfy the selection rules 



a + a 




1 mod 4 
3 mod 4 



l(a-a) E Z 



(7.5) 



Using the g-binomial building blocks and empirical selection rules, we thus obtain the fini- 
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Figure 13: Z 4 sectors (N, I odd): Combinatorial enumeration by single-columns of the q- 
binomial [^] = [g] = q~ 3 ^ 2 ^2 q^ 1 !m ' Ej . The excess of blue (even j) over red (odd j) 
1-strings is given by the quantum number a = \n/2\ — m = —2. The elementary excitation 
energy of a 1-string at position j is Ej = | (J — |) . The lowest energy configuration has energy 
E(a) = 1/4 + 5/4 = 3/2 = \o{a + |). At each empty position j, there is a 2-string. This 
analyticity strip is in the upper-half complex w-plane rotated by 180 degrees so that position 
j = 1 (furthest from the real axis) is at the bottom. The elementary excitations (of energy 
1) are generated by either inserting two 1-strings at positions j = 1 and j = 2 or promoting a 
1-string at position j to position j + 2. Notice that [™] ^ = [ n ^ m ] as g-polynomials but they 
have different combinatorial interpretations because they have different quantum numbers 
a. In the lower half-plane, q is replaced with q and no rotation is required. In this example, 
I = 7 and the value a = —2 of the quantum number in the lower half-plane is related to 
o = — 2 in the upper half-plane by the selection rules o + a = —(£ + l)/2 and |(er — a) e Z. 
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tized partition functions 

(tc)- 0/2J E«' 



Zf '(9) = < 



>2fc+£/2 



JV+1 
2 



- fc 





V-l 


g^2k-t/2 

- 1 


2 





fcez 

(gg)- c/24 E? A 



2fc+£/2 



2fc+£/2 



2fc+</2 



fcez 



JV+1 
2 

JV+|+2 + £ 

JV+1 
2 

^+2_£ 



JV+1 
2 



-A 



2k-C/2 



JV-1 
2 

N+t = 2_ k 



JV-1 
2 



iV-£ = mod 4 (7.6) 





JV-1 


q^2k-t/2 


2 


JV+l _ 


- <? 





N-£ = 2modA (7.7) 



For given mod 4 parities of N — £, these expressions are equivalent as partition functions 
but, in each case, the first form is used for the combinatorial interpretation when i = 1 mod 
4 and the second form when i = 3 mod 4. We further observe that 



£<iV 



JV+1 

r 2 



£e2N-l 



JV-1 JV+1 
2 2 
2n-l . -I— r , 2n-l , 



JV-1 
2 

2n — 1 % T — r , 2n— 1 . 



n=l 



n=l 



n=l 



n=l 



(7i 



7.2 Ramond sectors (TV even, £/2 even) 

The building blocks of the spectra in the Ramond sectors consist of the g-binomials 

--q- 1 ^ 2 E (f^ m]E \ a=[n/2\-m 

(7-doublc 
columns 



n 




n 


m 




\n/2\ - o_ 



(7.9) 



with Ej = j — ~. The sum is over all double-column diagrams as in Figure [15] with a fixed 
a. Here a is a quantum number given by the number of 1-strings in the right column minus 
the number of 1-strings in the left column 



a — m r ; g ht — mieft 

The number of 1-strings plus 2-strings at any given position is exactly two 

irtj + rij = 2, j = 1, 2, . . . , n 



(7.10) 



(7.11) 



The double-columns with quantum number a are generated combinatorially by starting 
with the minimum energy configuration of 1-strings for given a as shown in Figure [TU 
Empirically determined selection rules dictate that in a sector with I defects the quantum 
numbers of the groundstate satisfy 



a = a 



0,4,8, 



1 02 



(7.12) 



The energy of these groundstates is E(a) + E{a) — |(i . 

Excitations, incrementing the energy by one unit, are generated either by inserting a 
left-right pair of 1-strings at position j = 1 or incrementing the position j of a 1-string 
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Figure 14: Ramond sectors (N even, 1/2 even): Minimal configurations of double-columns 
with energy E(cr) = \o 2 . The quantum number a = \n/2\ — m is given by the excess of 
blue (right) over red (left) 1-strings. At each position j, the number of 1-strings rrij plus 
the number of 2-strings rij is 2. This analyticity strip is in the upper-half complex w-plane 
rotated by 180 degrees so that position j = 1 (furthest from the real axis) is at the bottom. 
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Figure 15: Ramond sectors (N even, 1/2 even): Combinatorial enumeration by double- 
columns of the g-binomial [^] = [g] = g -1 / 2 ^ qT,j m i E i m The excess of blue (right) over 
red (left) 1-strings is given by the quantum number a = \n/2\ — m = 1. The elementary 
excitation energy of a 1-string at position j is Ej = j — |. The lowest energy configuration 
has energy E(a) = \a 2 = |. At each position j, there are rrij 1-strings and rij = 2 — rrij 2- 
strings. This analyticity strip is in the upper-half complex w-plane rotated by 180 degrees so 
that position j = 1 (furthest from the real axis) is at the bottom. The elementary excitations 
(of energy 1) are generated by either inserting a left-right pair of 1-strings at position j = 1 
or promoting a 1-string at position j to position j + 1. Notice that \^\ q — [ n " m ] ? as Q~ 
polynomials but they have different combinatorial interpretations because they have different 
quantum numbers a. In the lower half-plane, q is replaced with q and no rotation is required. 
The value a of the quantum number in the lower half-plane is related to a in the upper half- 
plane by the selection rules a + a = 1/2 and \(o — a) G Z. 
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by 1 unit. The g-binomials are illustrated in Figure [15] Empirically, we find that all the 
excitations satisfy the selection rules 



a + a = £/2, 



la — a 



e z 



(7.13) 



Using the g-binomial building blocks and empirical selection rules, we thus obtain the frni- 
tized partition functions 



We further observe that 

KN 



■2k+e/2 



2[ N+2 ] 



£G4N 



q^2k-e/2 



2LfJ 



N even, f even (7.14) 



L n=l n=l 



n=l 



n=l 



(7.15) 



7.3 Neveu-Schwarz sectors (N even, i/2 odd) 

The building blocks of the spectra in the Neveu-Schwarz sectors consist of the g-binomials 



n 




n 


m 




\n/2\ - a 



q -\a{a+i) q ^i m i E i, a=[n/2\-m (7.16) 



cr-double 
columns 



with Ej = j. The sum is over all double-column diagrams as in Figure [T71 with fixed a. In 
these sectors, the number of 1-strings in the right column minus the number of 1-strings in 
the left column is related to the quantum number o by 



bright - "iieft = cr or a + 1 
The number of 1-strings plus 2-strings at any given position is exactly two 

rrij + rij = 2, j = l,2,...,n 



(7.17) 



(7.18) 



The double-columns with quantum number a are generated combinatorially by starting 
with the minimum energy configuration of 1-strings for given a as shown in Figure [16j For 
these minimum energy configurations 



bright — f^left — 0"min 



a, a>0 
a + 1, a < 



(7.19) 



Empirically determined selection rules dictate that in a sector with £ defects the quantum 
numbers of the groundstate satisfy 



a = a = {I -2)1 A, £ = 2,6,10,... 



(7.20) 
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Figure 16: Neveu-Schwarz sectors (N even, E/2 odd): Minimal configurations of double- 
columns within the binomials [^] = [ m ] ? - The energy is E(a) = \<j{o + 1) where the 
quantum number is o = \n/2\ — m. The excess of blue (right) over red (left) 1-strings in 
these minimal configurations is cr min as given in (17.191) . At each position j, the number of 
1-strings rrij plus the number of 2-strings rij is 2. This analyticity strip is in the upper-half 



complex u-plane rotated by 180 degrees so that position j 
is at the bottom. 
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Figure 17: Neveu-Schwarz sectors (N even, £/2 odd): Combinatorial enumeration by double- 
columns of the g-binomial [^] = [£] = q~ l ^ q^i m i E j . The number of positions is 
{n — l)/2 = 3. The quantum number is a = \_n/2\ — m = 1. The excess of blue (right) 
over red (left) 1-strings is a = 1 or a + 1 = 2. The elementary excitation energy of a 1-string 
at position j is Ej = j. The lowest energy configuration has energy E(a) = \cr(o + 1) = 1. 
At each position j, there are rrij 1-strings and rij = 2 — rrij 2-strings. This analyticity strip 
is in the upper-half complex u-plane rotated by 180 degrees so that position j = 1 (furthest 
from the real axis) is at the bottom. The elementary excitations (of energy 1) are gener- 
ated by either inserting a left or right 1-string at position j = 1 or promoting a 1-string 
at position j to position j + 1. Notice that ["] = [„ " m ] as g-polynomials but they have 
different combinatorial interpretations because they have different quantum numbers a and 
a' = —a — 1. In calculating a' , we have used the fact that n in ( 17.22ft is odd. In the lower 
half-plane, q is replaced with q and no rotation is required. The value a of the quantum 
number in the lower half-plane is related to a in the upper half-plane by the selection rules 
a + a = (£-2)/2 and \{a - a) G Z. 
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The energy of these groundstates is E(a) + E(a) = jg(^ 2 — 4). 

Excitations, incrementing the energy by one unit, are generated either by inserting a left 
or right 1-string at position j — 1 or incrementing the position j of a 1-string by 1 unit. 
The g-binomials are illustrated in Figure [TTJ Empirically, we find that all the excitations are 
consistent with the selection rules 



a + o 



(*-2)/2, 



a) ez 



(7.21) 



Using the g-binomial building blocks and empirical selection rules, we thus obtain the fini- 
tized partition functions 



^f } (g) = to)- c/24 E 



2Lf J + 1 ' 



2[ N±2 l _ 1 



J g 



N even, | odd 



(7.22) 



We further observe that 

KN 



LfJ 



e z i N) ^ = (^)"^n( i+ ^) 2 ri( i+ ^ 



n\2 



(7.23) 



n=l 



n=l 



7.4 Finitized characters in Ramond and Neveu-Schwarz sectors 



In the Ramond and Neveu-Schwarz sectors with N even, the (/-binomials do not give the 
appropriate characters. In fact, the relevant finitized irreducible characters are given [HI [25] 
by the generalized g-Catalan numbers 



n — 2r 
2 



chJ5(g) =q « E 



n-2 
2 



m=0 



m, m+r—l / q 



r+A 



1 -q r 

1 - q n l 2 



Here 

M \ 
m,n/ g 



n — 2r + l 
2 

! £«- 

m=0 



n-l 
2 



77, 

2 Jg 



-,2r 



c i A „ 1 — Z 
= (7-24+^2 - 

m,m+r — l/ q 1 — q n+1 



n even 



n + 1 

n+l 



L 2 



— r 



(7.24) 
n odd (7.25) 



lm(m+l)+|n(n+l) 



M 



M 
n 



- 9 



n— m+1 



M 
n + l 



M 
m — 1 



(7.26) 



are generalized versions [8] of the g-Narayana numbers [26] and the generalized Catalan 
numbers are 



.(») 

V,i 
.(») 



n-2 



r + iy v^-r-i 



ch^(l) 



n — 1 



n-l 



r + 1 



n — 1 



n-l 



r — 1 



r > 1, n even 



r > 1, n odd 



(7.27) 
(7.28) 



Despite the minus sign, (I7.26P is actually a polynomial with only non-negative coefficients. 
Combinatorially, the generalized g-Narayana numbers are generated by double-column dia- 
grams subject to admissibility (dominance) as shown in Figure [181 
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3 = 7 
j = 6 
j = 5 
j = 4 
3=3 
3 = 2 
3 = 1 



Figure 18: An admissible double-column diagram with m = 4 and n = 5 satisfying dom- 
inance L = {6,5,3,1} ^ {7,5,4,3,2} = R. The energy is E = 15 + 21 = 36 and the 
associated monomial is q 36 . 

A double- column configuration of height M consists of a pair of single-column configu- 
rations of height M. Suppose that there are m occupied heights in the left column and n 
occupied heights in the right column with < m < n < M. Let L = {L 1; L 2 , . . . , L m } be 
the occupied heights j in descending order of the left column and R = {R\, R2, • • • , R n } be 
the occupied heights j in descending order of the right column. We say the double-column 
diagram is admissible or satisfies dominance if L ^ R with the partial order 

L < R <^ L k < R k , j = l,2,...,m (7.29) 

Geometrically, one draws line segments between the occupied sites of greatest height in the 
two columns, then between the occupied sites of second-to-greatest height, and so on. Ad- 
missible double-column configurations are now characterized by not involving line segments 
with a strictly negative slope. 

We associate the monomial q E to each double-column configuration where 

m n 

e = j2 L k + J2 Rk ( 7 - 3 °) 

k=l k=l 

is the energy or weight of a double-column configuration. The generalized g-Narayana num- 
ber 

admissible 
double-columns 

is the sum of the monomials associated to the admissible double-column diagrams of height 
M and occupancy m and n in the left and right columns, respectively. 
The g-binomials decompose in terms of finitized irreducible characters 

n+l 

= £< +2) (?) (7-32) 

Q t=r 

n or n+l 

= E ch S n+1) (?) ( 7 - 33 ) 

1 t=r, by 2 




2n + 1 
n — r — 1 



r+A r o 



2n 
n — r + 1 
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Substituting the decompositions (I7.32p and (IT.33j) into (17.141) and (17.221) . we obtain the 
finitized partition functions 



E E Z >,'.' *S^ J+1, («) < 1 * J+1, (5), I even 



r=l r=l 



'(g) <K 



(7.34) 



E E^-< LtJ ^)< L ^"®' 



odd 



r=l r=l 

where the coefficients Z"^ r f are given by 



J £,r,r 



~(1 + (— l)' r+r ) [max (|, r + r) — max ( 



max (|, r + r) — max f ~, |r — r|), 



~, |r — r|j], | even 



odd 



(7.35) 



8 Conformal Partition Functions 



As N — > oo, the finitized characters become conformal characters. From (17. 6p . (17. 7p and 
(I7.34p . the conformal partition functions are given by the sesquilinear forms 



Ze(q) = < 



E ch 2 k+£/2(q)ch 2 k-£/2(q), ^4 {t odd) 

k=— oo 
oo 

^, r ,f ch r>2 (g) chr j2 (g), Ramond (| even) 

i-,f=i 

oo 

^, r ,f ch r i(g) ch fj i(g), Neveu-Schwarz (| odd) 



r,r=l 



Summing over sectors with either £/2 even or £/2 odd, we find 

'Zr,, = Mi + (-ir +p l [(t 1 ) + a 1 ) - I — 

I odd 



E^ 



(8.2) 



r r, 



and hence 



Z(q) 



E E cn 2fc+^/2(g)ch 2A; _£ /2 (g), Z 4 (£ odd) 

'G2N-1 fc=-oo 
oo 

Z rj f ch r>2 (g) ch f)2 (g), Ramond (| even) 

",f=i 

oo 

rf ch r) i(g) ch f) i(g), Neveu-Schwarz (| odd) 



(8.3) 



r,r=l 



Let us define the theta functions 



E * 



p\ 2 



(8.4) 
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It follows from (I7.6P and (17.71) . after summing over the Z4 sectors, that 

„ \#i 2 (<])\ 2 + \t)i 2 (<])\- 

£ Zt(q) = " m\>" (8 ' 5) 

where the Dedekind eta function is 

V (q) = q 1/2 *H(l-q n ) (8.6) 

71 = 1 

Similarly, from (I7.14p and (I7.22p . we note that summing over the I even sectors with suitable 
mutiplicities yields 

Z {q) + 2 2^ Zi{q) = -— = \x-i/s(q)\ + \XWW\ 

2 ^_ 2 Uq) ~ " W " 2|xo(g) + Xliq)l ^ 

1 3 

Z (q) + 2 J2 ZM) = rr-yp E l^(?)| 2 = \X-i/s{q)\ 2 + |X3/s(g)| 2 + 2|xo(g) + | 2 

where the W-irreducible characters are 

X-ys(q) = ^oM, Uq) = ^(qj^M+V(q) 3 } 

X3/s{q) = ^y^ 2 , 2 (g), xi{q) = frjjjftiM -v(q) 3 } 

In the last equation of (18. 7p . we have produced a modular invariant which we believe is 
the correct modular invariant partition function for critical dense polymers on the torus. 
However, we stress that this partition function is not obtained by the naive trace summing 
over all even N sectors since we have changed the multiplicity of all of the i 7^ sectors by 
a factor of 2. We note that the naive trace does not yield a modular invariant. 

Although we do not give details, the methods of this paper also apply to the case of 
Identified Connectivities (IC) so we include the results for completeness. In this case there 
are no defects, Ej = j and the finitized partition function is 

<'M = E<* J, M< ¥J) («> (8-9) 

r=l 

It follows that the conformal partition function is 

00 

z IC ( q ) = J2\cKi(q)\ 2 (8.10) 

r=l 

Further results could be obtained for the case of Distinguished Connectivities (DC) with £ 
defects and a = or a = —2. The cases ±2 are also integrable but can not be solved 

in the I = sector by the simple factorization techniques used in this paper. 
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9 Conclusion 



In this paper, we have solved a model of critical dense polymers exactly on arbitrary finite-size 
cylinders. This topology allows for non-contractible loops with fugacity a that wind around 
the cylinder or for an arbitrary number i of defects that propagate along the full length of 
the cylinder. We have set up commuting transfer matrices that satisfy a functional equation 
in the form of an inversion identity. For even N with loop fugacity a = 2, this involves 
a non-diagonalizable braid operator J and an involution R = — (J 3 — 12«7)/16 = (— 1) F 
with eigenvalues R = (— l)^ 2 . This is reminiscent of supersymmetry with a pair of defects 
interpreted as a fermion and F as a fermi number operator. However, we have no direct 
evidence of fermions or supersymmetry in our model. Nevertheless, since the number of 
defects £ is a quantum number that separates the theory into various sectors, we call these 
sectors Ramond {1/2 even), Neveu-Schwarz {1/2 odd) and Z 4 {£ odd), respectively. We 
conjecture that, when acting on link states with an arbitrary even number of defects, the 
braid operator R and transfer matrices T{u) are of rank 2 and have no rank-3 or higher 
Jordan blocks for arbitrary even N. The existence of non-trivial Jordan blocks is indicative 
that the theory on a long cylinder or torus is logarithmic in the sense that it gives rise to 
indecomposable representations of the Virasoro generator Lq. 

For a = 2 and arbitrary N, the inversion identity is solved exactly sector by sector 
for the eigenvalues of the transfer matrices in finite geometry. We thus obtain finitized 
conformal partition functions as sesquilinear forms in appropriate finitized characters. The 
eigenvalues are classified by the physical combinatorics of the patterns of zeros in the complex 
spectral-parameter plane. This yields selection rules for the physically relevant solutions to 
the inversion identity. The finite-size corrections follow from Euler-Maclaurin formulas. In 
the scaling limit, we obtain the conformal partition functions as sesquilinear forms and 
confirm [21 [3j H] the central charge c = —2 and conformal weights A, A = A t = {t 2 — l)/8. 
Here t = 1/2 and t = 2r — s G N in the I even sectors with Kac labels r = 1, 2, 3, . . . ; s — 1, 2 
while t e Z — | in the £ odd sectors. Strikingly, the £/2 odd sectors exhibit a W-extended 
symmetry but the £/2 even sectors do not. Moreover, the naive trace summing over all even 
N sectors does not yield a modular invariant. 

Recently, we considered [H] critical dense polymers on the strip with boundary condi- 
tions that, in the continuum scaling limit, respect a W-extended symmetry. Our analysis 
enabled us to identify this model with the c\p triplet model or symplectic fermions on the 
strip. It appears that no such simple identification can be made between these models on a 
long cylinder or torus after taking a naive trace. To resolve this seeming paradox, it may be 
necessary to modify the lattice model on a cylinder to properly restore the fermion structure 
and supersymmetry, or to introduce a modified trace implementing the proper torus geome- 
try. Only then can we hope to have a proper interpretation of the fermi number operator F 
and find agreement, in the N even sectors, with the conformal partition functions obtained 
by Saleur [I] . We intend to return to this question in a later publication. 
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A Properties of J 

For N even, J contains a non-contractible loop in exactly two diagrams 

J = J + Ji 

where Jq refers to the part without non-contractible loops, while 



i 

































(A.l) 



(A.2) 



(here illustrated for N — 8). In the enlarged TL algebra, this is just afl(E + Fj. 

In order to describe J, and not just J±, in terms of the enlarged TL algebra, we introduce 
the 'periodically ordered subsets' a n defined as follows. For n < N, let {01, . . . , a n } be a set 
of distinct elements of {0, . . . , N — 1}, and let a n be an ordering of {<Ti, . . . , a n } for which 
Oj+i 7^ &j — 1 and o"i ^ a n — 1. Two such orderings are equivalent if the corresponding 
ordered products of periodic TL generators match 



cr n ~ a n e ai . . . e c 



(A.3) 



We can now write J in terms of the enlarged TL algebra as 



N-l 



■i)^+E(-irEK-e,o) 2 + fi- 2 

n=l [ff„] 
2V-1 

n=l [ct„] 
2V-1 

•i)^ 2 + E(-^ n E K ■ ■ • e ^ -1 ) 2 + fi 

n=l [ct„] 
AT-1 

■i)^+E(- 1 ^E( rle -- e J 2+fi 



(A.4) 



n=l 



where the sum * s over the se ^ °f equivalence classes of periodically ordered subsets of 

cardinality n. a n is thus a representative of the given equivalence class. 



A.l Proof of Inversion Identity 

Proof of Inversion Identity (Section 14. ip Recalling that A 
corresponds diagrammatically to 



| , the left side of f l4T5|) 



7T , 



T(u)T(u + -) 



(A.5) 



with the left and right (vertical) edges identified. First, we examine the consequences of 
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having a half-arc connecting the two left nodes (or two right nodes) of a 3-tangle (two- 
column) appearing in (1A.5|) . Expanding in terms of connectivities, we find 



sin u cos u 



(^-f- +cos u sin u — sin 2 w + (j^- = — sin 2 w 



(A.6) 



IT) 



sin m cos u 



^ + cos it sin it ^) + -£) + cos 2 w ^~)= cos 2 w 



(A.7) 



and observe that such a half-arc will 'propagate' and ultimately lead to a vertical identity di- 
agram I. This is possible for a right- or left-moving half-arc. Taking the coefficients — sin 2 u 
(in the right-moving scenario (1A.6|) ) and cos 2 it (in the left-moving scenario (IA.7j) ) into con- 
sideration, this immediately gives the term proportional to J in (14. 5p . Having accounted for 
all situations with a half-arc, we are left with the two weighted connectivities 



— cos u sin u 



and 



cos it sin u 



(A.8) 



for each of the A" 3-tangles in ( 1A.5I) . Adding up all possible combinations of these connec- 
tivities readily produces the term proportional to J in (I4.5p . □ 



A. 2 Proof of Drop-Down Lemma 



Proof of Drop-Down Lemma (Section 14. 21) Since a non-trivial arc-part only appears for 
A^ > 2, we assume N > 2 and consider a segment of J consisting of two 3-tangles (14. ip with 
a half-arc linking the two nodes on the upper edge 





(A.9) 



The third diagram of the right side vanishes due to the contractible loop. For N = 2, the 
identification of the left and right vertical edges gives back a half- arc multiplied by 2 — a 2 , and 
the proof is complete. For A^ > 2, we include the 3-tangle to the left or right, respectively, 
of the first and fourth diagrams 






(A.10) 



(A.ll) 
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It follows that, for iV > 2, we simply have 



(A.12) 



implying that a single half-arc drops down leaving behind a couple of horizontal links (indi- 
cated by wavy links in (1A.12j) ) between the two vertical edges. For given arc-part of the input 
link state, this procedure is repeated starting with the innermost half-arc(s) and working 
out, as illustrated here 




in an example for N — 12. □ 

























































'4 














































(A.13) 



A. 3 Proof of Sector Lemma 



Proof of Sector Lemma (Section 14. 2p The lemma is readily verified for N = 1. For N = 2, 

there are two link states with £ = and one link state with £ = 2. «7 = Q 2 + Q~ 2 acts as 
twice the identity. J\ maps the link state with £ = 2 into minus the sum of the link states 
with £ = (which is outside the sector and therefore set to zero), while it maps both link 
states with £ = into themselves multiplied by —a 2 . This completes the proof for N = 2. 

N— £ 

For N > 2 and £ > 0, we get a factor of (— from the drop-down process. Within the 
given sector, the £ defects are mapped into £ defects. This requires the 3-tangles not affected 
by the drop-down process (as the 9th and 12th 3-tangles in the last diagram of (IA. 13[) ) to be 
identical thereby giving rise to an additional factor of ((—1)^ + 1). Since this factor is 2 for 
N, £ even and vanishes for N, £ odd, we immediately recognize ( 14. 7p . Finally, for N > 2 and 
£ = 0, we perform the drop-down process for all but the final, outmost half-arc. This gives 
a factor of (— l) - 2~~ and we are left with a scenario as in ( ]A.9j) giving rise to an additional 



factor of 2 



or 



□ 



A. 4 Matrix J for an arbitrary number of defects 

Here, we continue the description at the end of Section 14.21 of the matrix realization of J 
appearing in the Inversion Identity (1A.5I) for given iV when an arbitrary number of defects 
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is allowed. 

To simplify the description, let us introduce the following algebraic abbreviations 



u = -53 , d = E2 (A.14) 



allowing us to write a horizontal concatenation of such 3-tangles simply as products of it's 
and d's. For example, ignoring the half-arc linking the nodes on the upper edge in ( 1A.9I) . 
the sum of diagrams on the right side of flA.9[) reads uu + ud + du + dd. 



For a link state A, we let £(A) denote the number of defects, arc(A) the arc-part, and 
h(A) the maximal height of the arcs, where the height of an arc is the maximal number of 
points of intersection of arcs (on or below the arc itself) with a vertical line through the 
horizontal line segment between the two end-nodes of the arc. In the first diagram in ( 1A.13j) . 
there are thus three arcs of height 1, one arc of height 2, and one arc of height 3. 

Letting A and A' denote an input and an output link state, respectively, the corresponding 
entry of J is denoted by Ja',a- We immediately see that Ja>,a = if £(A) < £(A') since 
defects cannot be created by the action of J. Furthermore, the blocks with £(A) = £(A') 
appearing on the diagonal of J are themselves diagonal and given by the Sector Lemma 
(14. 7p . To complete the description of the matrix realization of J, we therefore turn to the 
cases where £(A) > £(A'). The Drop-Down Lemma implies that arc(A) C arc(A'), which we 
will assume in the following. As discussed in Section I4.2[ the data remaining to characterize 
the matrix J can be extracted by examining the input link state At, for which £(Ag) = I 
and arc(Ag) = 0, for each of the corresponding system sizes Ng = i where i = N mod 2. 

We let A' e denote an output state resulting when J acts on Ag. Since Ag consists of 
defects only, all arcs in A' e are created by J. This yields three cases distinguished by the 
height of A' e which can take on the values h(A' e ) =0,1,2 since Ja' ,A t = if h(A' e ) > 2. 

For h(A' e ) = 0, we see that A' e = Ag and J A ' v A e = (-1)' + 1- 

For h(A'g) = 1, we initially assume that £(A'g) = 0. This presupposes £ even, and we 
readily see that J r A' e ,A e = ( — l) 5 ^. For given A' e with £(A' e ) > 0, on the other hand, the 
contributing terms in the expansion of J have a specified distribution of 4-tangles ud 

. . . g%> {ud)g? . . . {ud)g? . . . g%> {ud)gf ... (A. 15) 

where every g is a u or a d, and where 

n 

£ = £' + 2n, J2 k i = £ ' ( A - 16 ) 

3=1 

Since £(A' e ) > 0, there is at least one segment between two neighbouring ucfs (which can be 
the same ud due to the cylinder topology) in which a u is not followed by a d, that is, n > 0. 
With reference to ( 1A.15I) . let us assume that this segment is the one indicated by g^ . . . g£\ 
If gjp = u, we see that g^ = u (g^ = u for n = 1) in order to prevent h(A' e ) = 2. In 

(2) 

order to preserve the specified distribution of (udys, every g\ in this segment must be a 
u. Continuing this argument around the cylinder, implies that all the g's must be it's. If, 
instead, g$ = d, this 3-tangle could be followed, on the other side of the (ud) to its right, 
by a g being a u. However, this u would force every g to its right to be a u, leading to an 
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inconsistency with g& = d once the cycle has been traversed. We thus conclude that all g's 

t(Aj)-t(A' t ) 

must be w's or all g's must be d's, implying that Ja' ,A t = (~ 1) 3 ((~ 1) + !)• 

For h(Af) = 2, there is at least one arc of height 2. In the contributing terms in the 
expansion of J, the 3-tangle containing the leftmost node of this arc must correspond to a 
u followed by a u to its immediate right, while the 3-tangle containing the rightmost node 
corresponds to a d preceded by a d to its immediate left, that is, 

. . . /£)(«« . . . dd)fU . . . fW( uu . . . dd)f[ 2) . . . fW(uu . . . dd)f[ 3) . . . (A.17) 

where every / is a u or a d, while Kj, v > 0. For v — 0, which presupposes I even, there 
is only one contributing term, and we see that Jma^ = ( — I) 3 - For v > 0, every segment 
ddfi . . . fniuu expands as 

<l<hlf ■ ■ ■ 9 { i X) {ud)gf 2) . . . glf . . . (ud)g[^ . . . g^fuu (A. 18) 

where ki, n > 0. As in the argument for h(A' e ) = 1, a g = u in flA.18j) cannot be followed by 
a g = d in f ]A.18j) . Ignoring the 4-tangles (ud), the sequence of g's in (1A.18j) can therefore 
be any of the Kj + 1 sequences of the form d . . . du . . . u. Adding up these contributions, 
while taking into account the 'boundary effects' coming from the sandwiching by dd . . .uu 
in (lA.18p . we see that all but one of the terms cancel out pairwise leaving 

1 + (— l) Kj 

dd(d . . . d(ud)d . . . d . . . (ud)d . . . djuu (A. 19) 

For v > 0, the 'value' for fj^) (that is, u or d) does not affect the value for fi +1 ^ (f[ for 
j — v — 1). This independence of the segments implies that J A' v A t is given by the product 
of the v contributions just obtained (being or 1) 

{0, if any Kj is odd 
1 -f 11 (A - 20) 
1, if ail Kj are even 

In summary, J = for odd, while, for N even, J is a sparse upper-triangular matrix 
whose non- vanishing entries are given by 

if A' = A 

if h A , A = 1 (A.21) 
if h A 'A = 2, s A 'a C 2N 



A', A 




where £' = £(A') and /ia'.a = h(arc(A') \ arc(A)), while sa>,a is the set of numbers of nodes 
between the rightmost node of a height-2 arc of arc(A')\arc(A) and the leftmost node of the 
next height-2 arc to its right (which, due to the cylinder topology, can be the same height-2 
arc). It is noted, that Ha\a = 1,2 requires that arc(A) C arc(A'). 

For even, we observe that the sum of the entries of a column of J is given by 



. „ a 2 -2, £ = 
C = -I T ^ (A.22) 
2a - 2, i = 2 mod N 



For given N, these sums are all equal if and only if a = 0, 2. 
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